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ASYMPTOTIC ZERO DISTRIBUTION FOR A CLASS OF 
MULTIPLE ORTHOGONAL POLYNOMIALS 

E. COUSSEMENT, J. COUSSEMENT, W. VAN ASSCHE 



Abstract. We establish the asymptotic zero distribution for polynomials gen- 
erated by a four-term recurrence relation with varying recurrence coefficients 
having a particular limiting behavior. The proof is based on ratio asymptotics 
for these polynomials. We can apply this result to three examples of multiple 
orthogonal polynomials, in particular Jacobi-Pineiro, Laguerre I and the ex- 
ample associated with Macdonald functions. We also discuss an application 
to Toeplitz matrices. 



1. Introduction 

Let /i be a positive measure on the real line for which the support is not finite 
and all the moments exist. The corresponding monic orthogonal polynomial P„ of 
degree n is then defined by 

(1.1) y"x™P„(a;)dAi(x) = 0, k = 0,...,n-l, 

with Pq = 1 and P_i = 0. A well-known fact is that such polynomials satisfy a 
three-term recurrence relation of the form 

(1.2) zPn{z)=Pn+l{z)+bnPn{z)+alPn-l{z), a„ > 0, 6„ € M, 

with initial conditions Po = 1 and P_i = 0. 

An object of frequent study is the asymptotic zero distribution of the zeros for 
a sequence of orthogonal polynomials. The zeros of the polynomials Pn, generated 
by (1.2), are real and simple [13]. With each polynomial P„ we can associate the 
normalized zero counting measure 

(1.3) y[Pn)--=- V <5., 

n ^-^ 

P„(x)=0 

where 5x is the Dirac point mass at x. If lim„^oo v{Pn) = v, by which we mean 
that 



lim [ fMPn) = f fdu 
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for every bounded and continuous function / on M {weak-k convergence) , then we 
call the probability measure u the asymptotic zero distribution of the sequence 
{Pn}n>o- One of the famous results in this context is the following. 

Theorem 1.1 (sec, e.g., [24, 28]). Suppose that the recurrence coefficients a„ > 
and 6„ e M have the limits a > and 6 e M, respectively. The polynomials Pn, 
generated by (1.2), then have the asymptotic zero distribution ^[7,5] with density 

n A\ d'^[7,<5]/ ^ / — \. a; e [7,(5], 

( 0, elsewhere, 
where 7 = 6 — 2a and S = b+2a. 

Remark 1.2. The measure W[-y,5] is known as the arcsine measure on [7, 6]. It also 
minimizes the logarithmic energy of the interval [7, S] [26]. 

Recently, the result in Theorem 1.1 was extended to the case of varying recur- 
rence coefficients. Here the notation lim„/jv^t Yn^N = Y denotes the property that 
in the doubly indexed sequence Yn,N we have lim^^oo Ynj,Nj = Y whenever Uj and 
Nj are two sequences of natural numbers such that Nj — > 00 and rij/Nj t as, 
j 00. 

Theorem 1.3 (Kuijlaars, Van Assche [20]). Let for each N G two sequences 
{on.Af}^!; On, TV > 0, and {fo„,Ar}J^Lo; bn^N G of recurrence coefficients be given. 
Furthermore, suppose there exist two continuous functions a : (0, +00) — > [0, +00), 
b : (0, +00) M, such that 

(1.5) lim an N = (i{t), lim bn n = b{t), t > 0, 

and define -f{t) b{t) — '2a{t), 5 := b{t) + 2a{t), t > 0. For the (orthogonal) 
polynomials generated by the recurrence 

(1.6) zPn,N{z) = Pn+l,Niz) + bn^PriMiz) + jv-^n-l.Af (-2^)) 

with initial conditions Po.n = 1 and P-i,jv = 0, we then have 

1 /■* 

(1.7) hm v{Pn,N) = - I ^[^{s),5(s)] ds, t > 0. 

n/N^t t Jq 

Here ^[7,5] is defined by (1-4) if'y<S and by 5-y if ^ = 5. 

Remark 1.4. More recently, Theorem 1.3 was generalized to measurable functions 
a and 6 [19]. 

In this paper we present a (conditional) theorem giving the asymptotic zero 
distribution for polynomials satisfying a four-term recurrence relation of the form 

(1.8) zPn,N{z) = Pn+l,N{z) + bn,NPn,N{z) + Cn,NPn-l,N{z) + dn,NPn-2,N{z), 

where the varying recurrence coefficients have some particular limiting behavior. So, 
in a sense it extends Theorem 1.3. Such a four-term recurrence relation appears in 
the theory of multiple orthogonal polynomials of Type II. These are a generalization 
of orthogonal polynomials which arises naturally in Hermite-Pade approximation of 
a system of (Markov) functions [9, 10, 22]. In particular, they satisfy orthogonality 
conditions with respect to several positive measures [2, 25, 29]. Some of their 
applications are situated in diophantine number theory, rational approximation, 
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spectral and scattering problems for higher-order difference equations and some 
associated dynamical systems, see, e.g., [5, 11, 21, 27]. Recently they also appeared 
in random matrix theory for matrix ensembles with external source [3, 7, 8] and 
Wishard ensembles [6]. The particular limiting behavior which we are considering 
appears in the examples Jacobi-Pineiro, Laguerre I [29] and the example associated 
with Macdonald functions [30]. 

In Subsection 2.1 we state our main theorem. Next, in Subsection 2.2 and Sub- 
section 2.3 we apply this result to the examples of multiple orthogonal polynomials 
mentioned above and some particular kind of Toeplitz matrices. In Section 3 we 
discuss a theorem on ratio asymptotics for monic polynomials satisfying the recur- 
rence (1.8). This will help us to prove our main theorem in Section 4. 



2. Statement of results 

2.1. Main theorem. We will study doubly indexed sequences of polynomials 
{Pn,N}, generated by a four- term recurrence of the form 

(2.1) zPn,N{z) = Pn+l,N{z) + b„^NPn.N{z) + C„^N Pn-l,N{z) + dn,N Pn-2,N {z) , 

with the initial conditions Pq.jv = 1, P-i,n = and P-2,n = and real recurrence 
coefficients. In particular, our main theorem gives an explicit expression for the 
asymptotic zero distributions 

lim i/(P„ jv), t>0, 

n/N^t 

with some conditions on the zeros of the Pn,N and some particular limiting behavior 
for the recurrence coefficients. As mentioned in the introduction the limit is taken 
over any sequence {v{Pnj,Nj)}j>i for which Uj oo, Nj oo and nj/Nj ^ x as 
j — > DO. We will use this notation throughout the rest of this paper. 

Theorem 2.1. Let for each N G N three sequences {bn,N}^=0' {cn,iv}^i and 

{dn.N}'i^=2 ^^al recurrence coefficients be given and assume that there exists a 
continuous function a : [0, +oo) — s- [0, +oo) such thai, for t > 0, 

(2.2) lim 6„,iv = 3/3(t), lim Cn,N=3Pitf, lim d„,Ar = /3(^)^ 

n/N^t n/N^t n/N^t 

with (3{t) = Let Pn,N be the monic polynomials generated by the recurrence 

(2.1) and suppose these polynomials Pn,N have real simple zeros x"'^ < . . . < x]^'^ 
satisfying the interlacing property x^~^^'^ < x^'^ < x^^^^'^, for all n,N € N, 
j = 1, . . . ,n. Then 

1 /•* 

(2.3) hm p{Pn,N) = - '"[o,a(s)] ds, t > 0, 

n/N^t t Jq 

where V[o^a] is defined by So ifa = and —^f^{x) = ^"~dx^(a)' ''^^^^ 



ifa>0. 



f V3 (i + VT^)V3 + (i-7r^)V3 

0, elsewhere, 
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0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



Figure 1. The density of the measure t'lo,!]- 



Remark 2.2. Denote by i^t the right hand side of (2.3). This measure acts on an 
arbitrary Borel set E like 

1 /•* 

iyt{E)^- V[„,^^,)]{E)ds, EeB{R). 
Jo 

Now suppose that for each a; > the set {s > | x < a{s)} is an interval, which 
we denote by [t-{x),t+{x)]. The density of the measure Pt is then 

,min(t,t+(x)) di;[o,a(.)] 
/ min(t,t_ (x)) 



(2.5) 



—(x)=- r 

dx ^""^ t L 



dx 



-{x) ds. 



This will be the case in each of the examples we present in this paper. 

Remark 2.3. Comparing Theorem 2.1 with Theorem 1.3 we see that the measure 
i;[o,i] plays the role of the arcsine measure in the case of orthogonal polynomi- 
als (satisfying a three-term recurrence relation). The density (2.4) again has the 
behavior ci(l — x)~^/^ as a; t 1, but has a different behavior C2a;~^/^ as a; J, 0. 

Remark 2.4. The measure W[o.i] coincides (after a cubic transformation) with the 
asymptotic zero distribution of Faber polynomials associated with the 3-cusped 
hypocycloid [18]. 



2.2. Application to multiple orthogonal polynomials. There are two types of 
multiple orthogonal polynomials but we will only consider type II. Let /xi, . . . , ^r, 
r G N, be a set of positive measures on the real line for which the support is not 
finite and all the moments exist. Furthermore, let n = (ni , n2, . . . , rir) be a vector of 

r nonnegative integers, which is a multi-index with length |n| := ni + n2 H h n^. 

A multiple orthogonal polynomial P^j of type II with respect to the multi-index 
n, is a (nontrivial) polynomial of degree < |n| which satisfies the orthogonality 
conditions 

(2.6) j x"'P.n{x)d^ij{x) =Q, < m < - 1, j = l,...,r. 

A basic requirement in the study of multiple orthogonal polynomials is that the 
system (2.6) has a unique solution (up to a scalar multiplicative constant) of degree 
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|n|. We call n a normal index for /ii, . . . if any solution of (2.6) has exactly 
degree \n\ (which implies uniqueness). If all the multi- indices are normal then the 
system of measures is called perfect. Some famous classes of perfect systems are 
the Angelesco systems, Nikishin systems ( for r = 2) and AT systems, see, e.g., 
[25, 29]. 

Multiple orthogonal polynomials of type II satisfy a recurrence relation of order 
r+1. In particular, if we set r = 2 and consider proper multi-indices Vn = {m+s, m), 
n G N U {0}, where n = 2m + s, s € {0, 1}, then the polynomials P„ := Pp^ satisfy 
a four-term recurrence relation of the form 

(2.7) zPn{z) = Pn+l{z) + bnPn{z) + CnPn-l{z) + dnPn-l{z), 

with the initial conditions Pq = 1, P_i s and P_2 = 0. For three examples 
known in the literature the recurrence coefficients in (2.7) have the particular lim- 
iting behavior (2.2), possibly after some re-scaling. In each of these examples the 
measures form an AT system on an interval A C R. It is then known that the 
zeros of the polynomials P„ are simple, lie in A [25, 29] and satisfy the interlacing 
property [4]. So, it is possible to apply Theorem 2.1. 

2.2.1. J acobi- Pineiro. The Jacobi-Pifieiro polynomials are the multiple orthogonal 
polynomials for the system of orthogonality measures 

dfj,j{x) = x"' (1 - xf dx, j = 1, 2, 

on the interval [0, 1] with ai, a2, /? > —1 and 0:2 — 0:1 ^ Z. In [29] it was shown that 

the monic Jacobi-Pificiro polynomials with respect to proper multi- indices, which 
we denote by P^^'"^''^, satisfy a recurrence relation of the form (2.7) for which 



(2.8) lim 6„ = 3 — , lim c„ = 3 — , hm rf„ = — 

By Theorem 2.1 with a{t) = 1, t > 0, wc then easily obtain the following result. 

Theorem 2.5. The Jacobi-Pineiro polynomials p^i'"^;/? /j^^g ^/jg asymptotic zero 
distribution f[o,i], defined as in (2.4). 

2.2.2. Multiple Laguerre I. The multiple Laguerre polynomials of the first kind are 
orthogonal with respect to the system of measures 

dfj,j{x) = a;"^e~'' da;, j = 1, 2, 

on [0, -l-oo) with oi, 02 > —1 and 02— Oi ^ Z. Denote the monic multiple Laguerre I 
polynomials with respect to proper multi- indices by L^^'"^. These satisfy a four- 
term recurrence relation of the form (2.7) where, for t > 0, 

bn 1 t \ , . Cn ^ 1 t \ , . dn ( t 



(2.9) lim = 3 - , lim = 3 - , lim 

see [29]. The following theorem is then a corollary of Theorem 2.1. 

Theorem 2.6. For the multiple Laguerre polynomials of the first kind the limit 

(2.10) u^:= lim V 6,/^, t > 0, 

L"1'"2(x)=0 
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exists and has the density 

(2-11) = 

"■^ 0, elsewhere, 

where, for y G (0, 1), 

3^/3 (1 + SVT^){1 - VT^)'^^ - (1 - 3^/^^)(l + ^/^^)^/^ 



giv) 



Proof. If wc define i"i'"^(z) := L"^ -""^ {N z) / N"^ , then the polynomials satisfy a 
recurrence relation of the form (2.1) and the asymptotic property (2.2) with a{t) = 
t > 0. So, applying Theorem 2.1 we get 

di/f 1 / 8 dujo,!] f8x\ _ 1 8 } I dt;[o,i] 



dx^ ^ t J 27s dx V 27s / 1 27 J u dx ^ ' 



Set a; € (0, ^) and ^; = Jl — Mi. Applying the substitution y ^ y/l — u we then 



27t 

obtain 



df,^, , 14\/3 

Hl + yy^/^l-yy^/^dy 

J —z 



-^{x) = i^y^ ((l + y)-^/^(l-2/)-^/^ + (l + y)-'^/^(l-y)-^/^)dy 

1 4^3 
~ t 277r 

1 73 l + 3y 



t IBtt (l+y)2/3(i_y)i/3 
This completes the proof. □ 
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2.2.3. Multiple orthogonal polynomials associated with Macdonald functions. In [30] 
one considered multiple orthogonal polynomials with respect to the orthogonality 
measures 

d/ii(.T) = x'^pj{x) dx, d/Lt2(x) = x'^p^+i(x) dx 
on (0, +00) with K > —1, 7 > and 

p^{x) = 2x'^/2i^^(2v/x), a; > 0, 

where is the modified Bessel of the second kind, also known as the Macdonald 
function [1, p. 374]. For the type II multiple polynomials with respect to proper 
multi- indices, P^''^, the recurrence coefficients in (2.7) are known [30, Theorem 4]. 
In particular, for i > 0, 

(2.12) hm 1^=3*2, lim = 3^^ lim = 

Theorem 2.1 then implies the following asymptotic result for the zeros of these 
polynomials. 

Theorem 2.7. For the multiple orthogonal polynomials associated with Macdonald 

functions the limit 

(2.13) := lim V 4/iv= 
exists and has the density 

, , / ^h{^), xe(o,^), 

[ 0, elsewhere, 

where 

Proof Define the polynomials PZ''^{z) := P2'"-{N'^ z)/N'^'' . By (2.12) these satisfy 
a recurrence relation of the form (2.1) having the asymptotic property (2.2) with 
a{t) = t > 0. Set a; e (0, and z = Applying Theorem 2.1 then 

gives 

du^' . 1 / A &v^( Ax\ 1 2 /• 1 

^("^ = 7 7 27;^^i27?j^^=i^277iy 7S^^"^^" 

2 /S" z 

3 V 3 

Similarly as in Theorem 2.6 we apply the substitution y <-> \/l — u and get 



1 1 1 + 



— Vl — 2 

Prom this we easily obtain (2.14). □ 
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2.3. Application to Toeplitz matrices. Set a > 0, (3 
Toeplitz matrices 



4a 
27 ' 



and define the 



/ 3(3 


1 





. \ 


3/32 


3/3 


1 






3/32 


3/3 1 ■• 









3/32 3/3 ■ • 


. 








1 



V 







n e N. 



/33 3/32 3/3 / 

Note that the eigenvalues of T" coincide with the zeros of the monic polynomial 
Qn{z) = det(2:/„ — T^), n G N. These polynomials satisfy the recurrence relation 

(2.15) zQ'^iiz) = + 3/3 + 3/32 Q«_i(^) +/33Q«_2(^), 

with Qq = 1 and Q"i = = 0- The following asymptotic result for the 

eigenvalues of the matrices T" then follows from Theorem 2.1. 

Theorem 2.8. The limiting eigenvalue distribution of the matrices T", with a > 0, 

is given by the measure ^[o.q], defined as in Theorem 2.1. 

Proof. The homogeneous recurrence relation 

= QJJ+i(0) + 3/3 0^(0) + 3/32 Ql_M + /3' Ql-2{^) 
with Q^(0) = 1 and Q^i(O) = Ql^iO) = 0, has the solution 

Q;:(0) = (-/3r(l + f + neN. 
So, since /3 = |y > all the are non-singular. Next, define 



/ 1 







3/3 


1 




3/32 


3/3 


1 




3/32 


3/3 







3/32 



\ 



^+l)x(n+l) 



: •. •• ■. •• ■•• 

V ... /33 3/32 3/3 1 y 

and notice that = {A'^f with 

/ 1 . 



A" = 



/3 1 
/3 



V 



■•• 

/3 1 / 



(n+l)x(n+l) 



Since is totally non- negative, T^ is totally non- negative, see, e.g., [15, p. 74, 
and so also T". By [15, p. 100, Theorem 10] we then get that the T" are oscillation 
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matrices. Consequently, the zeros of the polynomials Qn are simple and positive 
[15, p. 87, Theorem 6] and satisfy the interlacing property [15, p. 107, Theorem 14]. 
The theorem then easily follows from Theorem 2.1 and the recurrence (2.15). □ 

Remark 2.9. The polynomials satisfying the recurrence relation with constant 
coefficients (2.15) are the multiple Chebyshev polynomials of the second kind af- 
ter a cubic transformation [12, 14]. These are an example of multiple orthogonal 
polynomials of type II extending the well-known Chebyshev polynomials of the 
second kind [13]. The corresponding orthogonality measures can be found in [12, 
Corollary 4.2], [14, Theorem 4.1 ]. 



3. Ratio asymptotics 

In [20] Kuijlaars and Van Assche have proven a theorem that gives explicit ratio 
asymptotics for orthogonal polynomials with converging varying recurrence coeffi- 
cients. In this section we give an extension of this result to polynomials satisfying a 
(specific) four-term recurrence relation instead of a three-term recurrence relation. 

Theorem 3.1. Suppose we have for each N £N sequences {6„,jv}J^o, {cn,jv}$^i 
and {dn,N}^=2 of real recurrence coefficients and let Pn,N be the monic polynomials 
generated by the recurrence 

(3.1) zPn,N{z) = Pn+l,N{z) + bn,N Pn,N {z) + Cn,NPn-l,N{z) + dn,NPn-2,N{z), 

with Po,N ^ 1, P-i,N = and P-2,N = 0. Assume that for some fixed t > the 
recurrence coefficients have the limits 

\ / \ 2 / \ 3 

4a \ .. „ / 4a\ ,. , /4a 



(3.2) lim bn,N = 3 — , lim c„,jv = 3 — , lim dn,N , „„ 

n/N^t V 27 y n/N^t V 27 / n/N^t \ 27 

with a > 0. Furthermore, assume that the polynomials Pn,N have real simple zeros 
x"' < . . . < .t"'^ satisfying the interlacing property x^^^' < x"' < a;"^j^' , 
for all n,N e N, j = 1, . . . , n. Moreover, suppose that for some t* > t, there exist 
m < Q, M > a such that all zeros of Pn,N belong to [m, M] whenever n < t*N. 
Then 



(3.3) lim #^=<! 



i -*(-)• 

a Va/ 



a > 0, 



-, a = 0, 

z 



uniformly on compact subsets ofC\ [m, M] , where <f) is defined by 

(3.4) <p{z) := ^ {^-^ zV3 {^3(_i + 713^)1/3 + (_i _ ,/r37)i/3| _ 1 

with u)3 = e^ and 

\fpe^ = p^/^e^, p>0, 6le[0,27r). 

(3.5) (pe*^)i/3 = p'/^e^, p > 0, 0e(-^,+^]. 

Remark 3.2. In the case that the recurrence coefficients do not depend on A'^ the 
existence of the limit (3.3) was already proven in [4]. Our proof will be based on 
similar arguments. 
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Remark 3.3. Under the conditions of Theorem 3.1, by taking the derivative of 
(3.3) we get 



(3.6) 



1 ct>'{z/a) 



-Pr' + l,Ar(^) 



n/N^t\Pn,N{z) P„+i,Ar(z) 




a > 0, 
a = 0, 



uniformly on compact subsets of C \ [m, M] . 

In order to prove Theorem 3.1 we need part (a) of the following lemma. The 
whole lemma will be used in the proof of Theorem 2.1 as well. It can be found in, 
e.g., [20, Lemma 2.2] but we include a short proof for completeness. 

Lemma 3.4. Suppose that the zeros of the monic polynomials Pn-i and Pn, with 
degree n — 1 and n, respectively, are simple and real, interlace and lie in [m, M] . 

Then 



(a) 
(b) 



Pn-liz) 



Pn{z) 
Pn-l{z) 



< 



dist{z, [m,M])' 



yz eC\[m,M], 



Pn{z) 



- 2R' >max(|m|,|M|). 



Proof. Denote the real zeros of p„ by yi, . . . , t/„. Since Pn-i and p„ are monic and 
their zeros interlace, there exist Wj > 0, Y^^=i "^j = 1) so that 

Pn-l{z) 



Pniz) 



Then note that, because yj E [m,,M], for all z E C \ [m,M] wo have \z — yj\ > 
dist(z, [m, M]), 1 < j < n. This immediately proves part (a) of the lemma. 

If \z\ > max(|m|, |M|), then \yj/z\ < 1 and therefore 3?( j^_^^^ ) > |, 1 < j < n. 
Hence 



> 



Vj/z 

which proves part (b). 

We also need the following properties of the function (p. 
Lemma 3.5. The function (f> is analytic on C \ [0, 1] and satisfies 

mz) 



□ 



(a) zcPiz) = 1 + 



27 



^eC\[o,i], 



(b) = + 0(z"2), as z->oo. 

Proof. By the choice of branch cuts for the square and the cubic root, see (3.5), 
the function (j) is certainly analytic on C \ (— oo, 1]. For a; < 0, a simple calculation 
also shows that lim^j^o ^{x + ie) = lim^j^o <P{x — ie). So 4> is analytic on C \ [0, 1]. If 
we define 



(3.7) 
then 
(3.8) 



u±{z) :-- 
u+{z)u-{z) = z, 



u+{z) + u-{z) = —2. 
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Using this gives 



1 + 



40(£) 

27 



27 



27 



^(a;3(n+(^))V3 + (^_(^))i/3)' 



which verifies part (a) of the lemma. 

Next, we will prove that ^ tends to zero as z ^ oo. From part (a) we then obtain 
that lim^^oo zcpiz) = 1 and, since (p is analytic in a neighborhood of infinity, this 
implies part (b) of the lemma. Applying the formula (a + b) {a^ + 6^ — ab) = + b^ 
and (3.8), we observe that 

(3.9) ^6iz) = , 1. 

^7 u;| [(m+(z))1/3]2+ [(y_(^))l/3]2_^,^^l/3 

Now take for a moment z = 1 + L, with L > 0. By the definition of the square root 
we have ^±(1 + L) = —1 ± i\fL. So we can write 



w+(l + L) = p(L)e<^+^(^)), 
u_(l + L) = p(L)e^(-^-"(^)), 



< e(L) < 



Obviously we then get 



(3.10) 

(3.11) 
Finally, notice that 
(3.12) 



(«+(l + L))V3- 
(«_(1 + L))V3 



(p(L))2/3e'(-^+i^(^)), 
(p(L))2/3e'(-f-i^W). 



lim e{L) = and p{L) ~ \/Tj, L +oo. 



Hence, combining (3.9), (3.10) and (3.11), we obtain 

4 

lim — (^(1 + L) = lim 



-1 = 0. 



L^Too 27 '^'^ ' ' i^Tcx> gi(|£(L)) ^g-i(|£(L)) 

Since (I) is analj^ic in a neighborhood of infinity, then also lim^^oo (l>{z) =0. □ 

Now we give the proof of Theorem 3.1. 

Proof of Theorem 3.1: It is enough to prove the cases a = and a — 1. The 
more general case a > is then obtained by taking Pn,N{z) := Pn,N {az) /a". 

We first prove the case a = 1. By the assumptions on the zeros of the polynomials 
Pn,N every member of 



(3.13) 



\ Pn+l,N{z) 



n,N &n,n< t*N 



satisfies the estimate in part (a) of Lemma 3.4. So, the family (3.13) is uni- 
formly bounded on compact subsets of C \ [m, M] . By the theorem of Montel 
[17, p. 563] we then know that (3.13) is a normal family on C \ [m,M]. For a 
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sequence {(jij, Nj)}j>i, with nj,Nj oo, nj/Nj ^ t as j oo, we have that, if 
j is sufficiently large, the function 

(3.14) .m := J^"-^"-^l 

belongs to the normal family (3.13). The corresponding sequence {fj}j>i then has 
a subsequence that converges uniformly on compact subsets of C\ [to, M]. If we can 
prove that the limit of any such subsequence is (p, then, by a standard compactness 
argument, the full sequence {.fj}j>i converges uniformly on compact subsets of 
C \ [to, M] to (j). This then proves the theorem in the case a = 1. 

We will show that for each sequence ni,Ni oo with rii/Ni — » t such that the 
functions {/i}i>i converge uniformly on compact subsets of C \ [to, M], we have 

(3.15) f{z) := lim fi{z) = (j){z) + 0{z-^), as z ^ oo, 

for each fc G N. The uniqueness of the Laurent expansion around infinity then 
implies that f{z) = (p{z). We show this by induction on k. The case k = 1 follows 
from Lemma 3.5 (b) and fi{z) = 0{z^^), for every i > 1. Next, suppose that 
the claim holds for some k > 1 and consider a sequence {{rii, Ni)}i>i such that 
ni,Ni — > 00, rii/Ni x and the functions {/i}i>i converge uniformly on compact 
subsets of C \ [to, M] to some function / as z — > oo. If we put 

9i{z) ■■= zeC\[m,M], 

hi{z) := p-^-^'^'i^j , ^eC\[TO,M], 

then from the recurrence relation (3.1) we obtain 

(3.16) z = Mzy^ + bn^^Ni + CmM, gi{z) + rf«.,Ar. hi{z). 

Since t < t* we may assume without loss of generality that rii < t*Ni for every i > 1. 
Then {5fi}i>i and {/ii}i>i are subsets of the normal family (3.13). Therefore, there 
is a sequence ij oo, j oo, such that {gij}j>i and {hj.}j>i converge uniformly 
on compact subsets of C \ [m, M] with limit g and h, respectively. If we pass to 
such a subsequence and take limits in (3.16), then by (3.2) we find 

1 /4\ /4\^_ /4^^ 



z ^ oo. 



(3.17) ^ = _+3(^-J +3(^-J 9{z)+\^-] g{z)h{z), zeC\[m,M]. 

By the induction hypothesis we now have that 

g{z) = <^(^) + 0(^-'=), 

h{z) = (t>{z) + 0{z-''), z^oo. 
Applying this to (3.17), by Lemma 3.5 we then get 

4 \ o 4 , / 4 V 



= ^-3 - -3 - 0(z)-h^ cl>{zr + 0{z-'^) = — + 0{z-'^). 



f{z) \27 J \27 J ' V27y ' ' ' ' (t>{z) 

Since (/)(z) = 0{z~^), this implies 

^^'^ - l + ^{z)0{z->^) - l + 0(.-^-i) - -^("^ + 
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So we proved that (3.15) also holds with k replaced by A; + 2. Therefore, it holds 
for all k. 

Finally, for the case a = the proof is similar. In fact, (3.3) then easily follows 
by taking limits in (3.16). □ 

4. Proof of Theorem 2.1 

In order to prove the asymptotic result in Theorem 2.1 we first have a closer 
look at the function which is analytic on C\ [0, 1]. Here we will use the relation 



Hz) ^ l^s^ _ -^^ 2z 2z \ 27 



which can be obtained by differentiating part (a) of Lemma 3.5. First of all, we are 
interested in the jump across its branch cut. 

Lemma 4.1. The jump of the function ^ across its branch cut is given by 

(4.2) mix) := lim \) . ,^ - lim = 2m —f^ix), a;e(0,l), 

where f[o,i] is defined as in Theorem 2.1. 
Proof. Let x G (0, 1). By (4.1) we easily obtain 

(4.3) m(x) = — lim ' - 1 - lim -^-^ — ^ - 1 
^ ^ ^ ^ 2a; [eio V 27 J eio \ 27 

Applying the definitions (3.5) we get 

1 — 1 + V i — xe^'- ' ' ^ ^ 

eJ.0 



lim(-l + \/l-xe±^'^) = e't (l± Vl-a;)^ 
im (^-1 - V 1 - a;e±»^ j = e"* 3 (l ip Vl - x) 



eJ.0 

Using the notation v±{x) := (1 ± ^/l — x)^^^ we then have 



lim ■ 



;xe±-) 



1 = 3a;i/3 |e-¥i;^(a;) + e"? u=p(x)| - 3. 



eio 27 

So, also applying the relations x^^^ = v+{x)v-{x) and 2 = w+(a;)^ + V-{x)^, equa- 
tion (4.3) becomes 

2V^i v^[x) — V-[X) 

a;^/^ {x^/3{v+{x) + v^{x)) - 2f + 3,x2/3 {v+{x) - ^-(x))^ 

= '^-^{v+{x)-v.{x))-' [{v.{x)^ -v+{x)^)^ +^v+{x)^v.{x)^]-\ 

If we multiply the numerator and denominator both by v+(a;) + V-{x), then we 
finally obtain 

2\/2>i v+[x)+V-{x) v^i v+{x)+V-{x) 

= a;2/3 {v+{x)^ + v-{x)^) (i;+(x)3 - v^ix)"^) " 2^ ^T^^ ' 

which proves (4.2). □ 
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A second point of interest is the behavior of the function ^ in the neighborhood 
of its branch points and 1. 



Lemma 4.2. Near the points and 1 we have 



(4.4) 
(4.5) 



)(£e*^) 
/>'(!+ ee^^) 



= O 



o 



,-2/3 



elO, 



)(1 +ee*^) 

Proof. For 6 € (0, 27r) we easily see that 



0e (0,2^), 

6 e (-TTjTr) . 



27 

Um<^(£e*^) = -- 

ej.0 ^ ' 4 



4^(ee'^) 



27 



+ 1 



Applying this to the first equality in (4.1) then gives expression (4.4). 

We now take 9 € (— 7r,7r). We then have u±{l + ee*^) = — 1 ± y/se^^ , where 
we use the notation (3.7), and applying the definition of the third root (3.5) we get 

(4.6) hm</,(l +ee^«) = f {c^ge^ +e-^} - l) = |. 

If we write u±{l +ee'^) := p±{s,9)e^'^^^^'^\ meaning the polar coordinates, then a 
closer look gives 

ri±{e,0)=±TrTV£cos{9/2)+O{e), e i 0, 

and 



(4.7) 



.'7±(=,«) 



we 



p±ie, 9)'/^ = 1 ± ^ sin (9/2) + 0(e), e i 0. 



(4.8) 

From this we easily obtain 
(4.9) + 



27 



- 1 = -Vsie^ +©(£), £iO. 



Applying (4.6) and (4.9) to the first equality in (4.1) then finally leads to (4.5). □ 

As a corollary of Lemma 4.1 and Lemma 4.2 we obtain that ^ is the Stieltjes 
transform of the measure f[o,i], up to a minus sign. 

Lemma 4.3. Let be defined by (3.4), then 



(4.10) 



£(£) 



[ ^—dv[o,i]{x), zeC\[0,l], 
J z X 



with W[o,i] defined as in Theorem 2.1. 

Proof. By Lemma 4.1 and Lemma 4.2 and applying Lemma 3.5 (b) to (4.1), the 
function ^ satisfies the following additive Riemann-Hilbert problem: 

(PI) / is analytic in C \ [0, 1], 

(P2) lim/(a; + i£) -lim/(a;-ie) = 27rz ^^(x), for (0,1), 



e4.0 



£4.0 
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(P3) f{z) = -z-^+0 (z-2) , as z ^ oo, 



(P4) f{z) = O 

f{z) = 0{{l-z)-'/^) 



as z 0, 
as z 



If / and g are both solutions of this Ricmann-Hilbcrt problem, then it is easily seen 
that f — g is analytic in C\{0, 1}. Moreover, and 1 are removable singularities by 
(P4). Liouville's Theorem and (P3) then imply f = g, meaning that the Riemann- 
Hilbert problem has a unique solution. So it is enough to show that 

f{z) := I ^— di;[o,i](a;), ^ e C \ [0, 1], 
J X z 



satisfies (P1)-(P4). Properties (PI) and (P3) easily follow from the fact that 

is a probability measure on (0,1). By the Plemelj-Sokhotskii formula for Cauchy 

integrals, see, e.g., [23, p. 43, (18.1)], / satisfies (P2). Finally, the behavior at 
the branch points, see (P4), easily follows from [23, p. 74, (29.5) and (29.6)] and 
(2.4). □ 



Remeirk 4.4. As an easy consequence of Lemma 4.3 we obtain 

1 (l)'iz/a) 



(4.11) 



z — X 



a (j){zla) 
1 



a > 0, 
a = 0, 



z e C \ [0,a], 



where V[o,q] is defined as in Theorem 2.1. 



We are now ready to prove Theorem 2.1. 

Proof of Theorem 2.1: Let f > and fix a number t* > t. Clearly, the con- 
vergence (2.2) and the fact that the function a is continuous on [0, oo) imply that 
the recurrence coefficients are uniformly bounded if n/N is restricted to compact 
subsets of [0, oo). So, 

(4.12) Q<R:= sup{l + |6„,jv| + |c„,jv| + |rf„,jv| : n < t*N} < +00. 

By the recurrence (2.1) we have Pn,N{z) = det(^;/„ — Ln,N), with 



(4.13) 



N 



1 bo,N 


1 





\ 


Cl,N 


&1,JV 


1 




d2,N 


C2,N 


b2,N 1 







d3,N 


C3,N ^3, AT 





V 




dn-l,N 


1 

Cn-1,N bn-l,N / 



As a consequence, the zeros of Pn,N arc bounded by ||in|loo- For n < t*N we then 
have that the simple zeros of P„,jv lie in the interval [—R,R]. Moreover, they are 
assumed to satisfy the interlacing property (for fixed N) and one can observe that 



a{st) < 1 + 



4a(si) 
27 



< R, 



for each < s < 1. 
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So, by Remark 3.3 and Remark 4.4 we establish 



(4.14) 



lim 



P[sn\A^)\ / 1 ^ 
PYsn\,N{z) j J Z-X 



i(st)]{x) 



uniformly on compact subsets of C \ [—R,I{\, where < s < 1 and [snj denotes 
the greatest integer less than or equal to sn. 

Note that 
(4.15) 



1 P^(^) _ 1 



P'k+iiz) P',{z) 



nPn{z) n^^\Pk+l{z) Pk{z)J Jo yPlsni+l{z) P[sn\{z) 



I 

Jo 



P[ 



{z) 



P[sn^ 



ds. 



For n < t*N the zeros of the polynomials Pn,N arc simple, lie in [— i?, -R] and satisfy 
the interlacing property for fixed N . From Lemma 3.4 (b) we then get 

Pysn\+l,N{z) 



< %z\ 



\z\ > R. 



Plsn\,N{z) 

With a similar argument as in Lemma 3.4 (a) we can also prove 



Plsn\,N{z) 

,n{z) 



LsnJ + 1 



< 



dist(^,[-i?,i?])2 



zeC\[-R,R]. 



Combining these two results we have, for l^] > i?. 



LsnJ + 1 



{z) 



[snJ 



{z) 



Pysn\+l,N{z) 



P'lsn\,N{z) 

2\z\ 



,jv(^) 



P'lsn\+l,N{z) 



^^'^^^ - dist(z, hi?,i?])2' 

So, we can apply Lebesgue's dominated convergence theorem on (4.15) and by 
(4.14) we obtain 



lim 

n/N^t 



1 



■ A.v{Pn,N){x) 



1 



(4.17) 



z — X 
1 



z — X 



dw[o,a(*it)] (a;) ds 
du[o,a(s)](a;)ds. 



for 1^1 > R. By [16, Theorem 2], which is a gloss on the theorem of Grommer and 
Hamburger [31, p. 104-105], we then finally establish (2.3). □ 
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